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Abstract. The High Energy Stereoscopic System (H.E.S.S.) has observed the high-frequency peaked BL Lac object 
PKS 2155-304 in 2003 between October 19 and November 26 in Very High Energy (VHE) y-rays (E > 160 GeV for these 
observations). Observations were carried out simultaneously with the Proportional Counter Array (PCA) on board the Rossi 
X-ray Timing Explorer satellite (RXTE), the Robotic Optical Transient Search Experiment (ROTSE) and the Nancay decimet- 
ric radiotelescope (NRT). Intra-night variability is seen in the VHE band, the source being detected with a high significance on 
each night it was observed. Variability is also found in the X-ray and optical bands on kilosecond timescales, along with flux- 
dependent spectral changes in the X-rays. A transient X-ray event with a 1500 s timescale is detected, making this the fastest 
X-ray flare seen in this object. No correlation can be established between the X-ray and the y-ray fluxes, or any of the other 
wavebands, over the small range of observed variability. The average H.E.S.S. spectrum shows a very soft power law shape 
with a photon index of 3.37 ± 0.07 stat ± 0. 10 sys . The energy outputs in the 2-10 keV and in the VHE y-ray range are found to be 
similar, with the X-rays and the optical fluxes at a level comparable to some of the lowest historical measurements, indicating 
that PKS 2155-304 was in a low or quiescent state during the observations. Both a leptonic and a hadronic model are used 
to derive source parameters from these observations. These parameters are found to be sensitive to the model of Extragalactic 
Background Light (EBL) that attenuates the VHE signal at this source's redshift (z = 0.117). 

Key words. BL Lacertae objects: Individual: PKS 2155-304 — galaxies: active — gamma rays: observations — X-rays: 
galaxies — radiation mechanisms: non thermal 



1. Introduction 

The innermost regions of active galactic nuclei, where the 
largest part of their luminosity is emitted, can be probed 
through observations of their flux variability at different wave- 
lengths. The physical processes in their central engines and 
jets are usually considered the main candidates for the origin 
of the observed variability. Measurements of correlated vari- 
ability, spectral variations and time lags across the broad-band 
observations allow modelling of particle distributions and their 
radiation processes, as well as probing the acceleration mecha- 
nisms that are involved. 

PKS 2155-304 is probably the most prominent and best- 
studied blazar-type Active Galactic Nuclei (AGN) in the 
Southern Hemisphere. The emission of PKS 2155-304, its 
possible variability patterns, as well as correlations across all 
wavebands, have been studied exhaustively over the past 20 
years (see e.g. |Urry et al. 1997) . Its first detection at VHE y- 
rays by the Durham Mk VI telescopes ( Chadwi ck et al. 1 999 ) 
classified it as a TeV blazar, like the northern hemisphere 
BL Lac objects Mkn 421, Mkn 501, H 1426+428, or 
1ES 1959+650. Its redshift of z = 0.117 makes it the sec- 
ond most distant confirmed TeV blazar after H 1426+428 (z = 
0.129). PKS 2155-304 was the brightest BL Lac object in the 
EUVE all-sky survey (Marshall, Fruscio ne & Carone 1 995 1. 
This source was confirmed as a high energy y-ray emitter by 
H.E.S.S. JAharonian et al. 20041 AH04 hereafter) at the 45<x 
significance level, when strong detections were reported for 
each of the dark periods of observations. 

Here we report on simultaneous H.E.S.S. VHE y-ray, 
RXTE/PCA X-ray, ROTSE optical, and NRT decimetric 
observations of PKS 2155-304 during the dark periods 
of October and November 2003. No simultaneous multi- 
wavelength campaign had before included an Atmospheric 
Cherenkov Telescope (ACT) that could sample the evolution of 
the high energy component of the spectral energy distribution 
(SED) of this object. We also include EGRET archival data, 
and other archival radio through X-ray data obtained from the 
NASA/IPAC Extragalactic Database (NED). Details of the ob- 
servations and data reduction/analysis are given in ^2] Light 
curves and spectra are described in ^3] The attenuation of the 
H.E.S.S. spectrum by the EBL and an interpretation of the data 
using a leptonic and a hadronic model are discussed in 

2. Observations and data analysis 

2.1. H.E.S.S. 

2.1.1. H.E.S.S. detector 

In its first phase, the H.E.S.S. array consists of four atmo- 
spheric Cherenkov telescopes operating in stereoscopic mode. 
However, the data shown here were taken during the construc- 
tion of the system, initially with two telescopes, with one more 
as of mid-September, 2003. The fourth and final telescope was 

Send offprint requests to: Berrie Giebels e-mail: 
berrie@poly . in2p3 . fr 

* UMR 7164 (CNRS, Universite Paris VII, CEA, Observatoire de 
Paris) 



added to the array in December 2003, subsequent to the data 
presented here. Each telescope has a tessellated 13 m-diameter 
( 1 07 m 2 surface area) mirror which focuses the Cherenkov light 
from the showers of secondary particles created by the inter- 
action of y-rays in the atmosphere onto a camera in the focal 
plane. This camera consists of 960 photomultipliers, each with 
a pixel size of 0.16°, giving a field of view of 5°. For the data 
sets presented here the y-ray trigger threshold is m 100 GeV 
and the spectral threshold is 160 GeV with an energy resolu- 
tion =a 15%. The experiment is located in the Khomas high- 
lands in Namibia, (23° S, 15° E, 1800 m a.s.l.). A detailed de- 
scription of the layout and the components of the telescope op- 
tical systems, including the segmented mirror with its support 
structure, the mirror facets for each telescope and the Winston 
cone light concentrators in front of the PMT camera can be 
found in Ber nlohr et al. (2 03) | and a descript ion of the mir- 
ror alignment is given in IConul^^L^2003)|^F2or details on 
the camera calibration see Aharonian et al. (2004) The trigger 
system is described in |Funk, Hermann, Hinton et al. (2 004) 
Early reports of H.E.S.S. have been given elsewhere (see e.g. 
IHofmann 20021 . 

2.1 .2. Observations 

In the beginning of the observation period of October 2003, 
a single 5 cr detection by H.E.S.S. achieved in lh, triggered 
an approved RXTE ToO (Target of Opportunity) on this target. 
The results presented here are based on observations carried 
out in 2003 between October 19 and November 26. 

The data were taken in the Wobble mode where the source 
direction is positioned ±0.5° in declination relative to the cen- 
tre of the field of view of the camera during observations. This 
allows for both on-source observations and simultaneous esti- 
mation of the background induced by charged cosmic rays. The 
data reported here are selected and analyzed with the "standard 
analysis" described in §4 of AH04. The background is esti- 
mated here by using all events passing cuts in a ring around 
the source, as described in section 4.3 in AH04. The runs pass- 
ing the quality selection criteria total 32.4 hours of livetime on 
the source. The total two-dimensional significance sky map is 
shown in Figure^ along with a graph showing the 6 2 distribu- 
tion (where 6 2 is the square of the angular difference between 
the reconstructed shower position and the source position) of 
the 1764 excess events observed. This yields a detection at the 
34. 3<t level, at the average rate of 0.9 ly/ min and a significance 
of 6.0o7 Vh. 

The methods used here for reconstructing the energy of 
each event and for determining a spectrum are described in §6 
of AH04. The measured time-average spectrum is fitted by a 
power law of the form dN/dE = Io(E/ 1 TeV)~ r with Iq the flux 
normalization at 1 TeV and T the photon index. The photon in- 
dex obtained from the time-averaged spectrum is then used as a 
fixed parameter to estimate the integral flux above 300 GeV for 
each run. This integrated flux takes into account the effective 
area and threshold variations due to the source moving through 
the sky, giving more reliable variability information than count- 
ing rates in units of y-ray s/minute. Overall systematic errors 
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vided by NASA/GSFC and filtered using the RXTE Guest 
Observer Facility (GOF) recommended criteria. The changing 
Proportional Counter Units (PCUO/2/3) configuration through- 
out the observations was taken into account in the data reduc- 
tion, and only the signals from the top layer (X1L and X1R) 
were used. When reducing the PCUs individually to estab- 
lish the time-averaged spectrum, the average spectral fit pa- 
rameters are similar within error bars but with a systemati- 
cally higher^- 2 for spectral fits performed on PCU0 data alone. 
On shorter timescales this effect becomes negligible and PCU0 
contributes to the statistical significance of the flux measure- 
ment. Therefore all PCUs were kept for the overall light curve 
which is binned in 400 s bins, but only PCU2 and PCU3 are 
used for the analysis of data segments that are simultaneous 
with H.E.S.S. runs and for the time-averaged spectrum. 

The average spectrum used in the SED is derived by 
combining PCU2 and PCU3 spectra using the addspec tool 
weighted by the counts information delivered by f statistic 
and then the corresponding response matrices were combined 
with addrmf . The faint-background model was used and only 
the 3-40 PHA channel range was kept in XSPEC v . 11.3.1, 
or approximately 2-20 keV. 

To build a light curve in units of integrated flux in the 
2-10 keV band, spectral data were derived from 400 s bins 
to probe short timescales with adequate statistical accuracy. 
These segments are then fitted by a power law in XSPEC 
with PCU configuration-dependent response matrices gener- 
ated by the ftool pcarsp v. 10. 1 and a fixed column den- 
sity of N H = 1.7 x 10 20 cnT 2 obtained from PIMMS 1 . This 
yields the flux and the error (corresponding to the lcr confi- 
dence level) on the flux reported in the light curves in Figure|5J 
in units of 10~" erg cirT 2 s _1 in the 2-10 keV band. The fits did 
not improve by using a broken power-law for the 400 s binned 
observations. 
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Fig. 1. Top: The distribution of 6 2 for on-source events (points) 
and normalized off-source events (shaded) for the observations 
in October and November 2003. Bottom: The two-dimensional 
distribution of excess events observed in the direction of PKS 
2155-304. The bins are not correlated and represent the actual 
distribution of observed y-rays on the sky. The right hand scale 
is the number of counts. 



are estimated to be 20% for the integral flux and ^0.1 for the 
photon index. 



2.2. RXTE 



The PCA dJahoda et al. 1996i units of RXTE observed 
PKS 2 1 55 -304 between October 22 and November 23 of 2003 
with exposures of typically lOks in October and 1 ks in 
November. The STANDARD2 data were extracted using the 
ftools in the HEASOFT 5.3.1 analysis software package pro- 



2.3. ROTSE 

The ROTSE-III array is a worldwide network of four 0.45 m 
robotic, automated telescopes built for fast (« 6 s) response to 
GRB triggers from satellites such as HETE-2 (High Energy 
Transient Explorer 2) and Swift. The ROTSE-III telescopes 
have a wide ( 1 .85 x 1 .85 °) field of view imaged onto a Marconi 
2048 x 2048 pixel back-illuminated thinned CCD and are op- 
erated without filters. The ROTSE-III systems are described 
in detail in Akerlof et al. (2003) At the time of the observa- 
tions of PKS 2155-304 in October and November 2003, two 
ROTSE-III telescopes were operational in the Southern hemi- 
sphere: ROTSE-IIIa located at the Siding Spring Observatory, 
Australia and ROTSE-IIIc at the H.E.S.S. site. The ROTSE- 
IIIc telescope is located in the centre of the H.E.S.S. telescope 
array. A 30% share of the total observation time is available 
to the H.E.S.S. collaboration, which has been used to perform 
an automated monitoring programme of blazars, including ob- 
jects that are being observed with the H.E.S.S. telescopes. 
Both telescopes participated in the observation campaign on 
PKS 2155-304 in October and November 2003. 



See http://legacy.gsfc.nasa.gov/Tools/w3pimms.html 
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The telescopes observed PKS 2155-304 typically 10 times 
per night taking sequences of 2 frames with 60 s exposures with 
a slight dithering of the pointing to reduce the impact of indi- 
vidual noisy pixels. The typical limiting magnitude, depending 
on the sky conditions, is 18.5 mag . Overall, 323 bias-subtracted 
and fiat-fielded frames have passed visual inspection and are 
used to produce a light curve. A total of 6 frames were rejected 
due to the presence of stray light from Jupiter. 

Using an overlay of 50 isolated comparison stars with 
similar brightness (12 mag - 14 mag ) and co-located with 
PKS 2155-304 (< 15 arcmin), a two-dimensional Gaussian is 
fit to the intensity distribution characterising the point-spread 
function by cr ps f . 

To estimate the local sky-background for the reference stars 
and the target object, an annulus with inner radius 2 x cr ps f and 
an outer radius 6xcr ps f is chosen. Based upon a reference frame 
which is derived from co-adding 30 individual frames, a mask 
is calculated for each object excluding regions where faint ob- 
jects coincide with the annulus; pixels exceeding 3 standard 
deviations of the local sky background are excluded. Using the 
local sky background, the intensity and error of each object is 
calculated. Using the 50 reference stars, a relative intensity and 
statistical error with respect to a reference frame is calculated. 

The absolute flux values are obtained by calculating a rel- 
ative R magnitude by comparing the instrumental magnitude 
with the USNO catalogue as described in |Akerlof et al. (20 00) 
The procedure has been checked by comparing the average R 
magnitude of a sample of 70 BL Lac type objects determined 
with ROTSE observations carried out over one year of oper- 
ation with the V magnitude listed in the 10 th Veron Cetty & 
Veron catalogue of BL Lac type objects. The average V - R of 
0.5 that is found is consistent with the average value for V — R 
obtained from cross-checking the colours with the 2MASS cat- 
alogued value for the BL Lac type objects. 

Finally, the host galaxy has been resolved in optical 
JFalomo 1996l > and NIR ( |Kotilainen, Falomo, & Scarpa 1998| 
KFS98 hereafter) and found to be a typical giant elliptical with 
M(R) = -24.4 which translates into an apparent m(R) = 15.1 
(here the distance moduli given by KFS98 have been used to 
calculate the apparent magnitude based upon the absolute mag- 
nitude quoted). The ROTSE measurements have as maximum 
and minimum m(R, min) = 13.3 and m(R, max) = 13.7 which 
corresponds to lOmJy for the maximum observed flux and 
6.7 mJy for the minimum flux taking the contribution of the 
host galaxy into account. These values are considerably lower 
than the retrieved archival data indicating that PKS 2155-304 
was in a low state at the moment of the observations. 

2.4. NRT 

The Nan§ay radiotelescope is a single-dish antenna with a col- 
lecting area of 200 x 34.56 m 2 equivalent to that of a 94 m- 
diameter parabolic dish (van Drie l et al. 199 61. The half -power 
beam width at 1 1 cm is 1.9 arcmin (EW) xl 1.5 arcmin (NS) (at 
zero declination), and the system temperature is about 45 K in 
both horizontal and vertical polarizations. The point source ef- 
ficiency is 0.8 K Jy~', and the chosen filter bandwidth was 12.5 



MHz for each polarization, split into two sub-bands of 6.25 
MHz each. Between 4 and 14 individual 1 -minute drift scans 
were performed for each observation, and the flux was cali- 
brated using a calibrated noise diode emission for each drift 
scan. Data processing has been done with the Nancay local 
software NAPS and SIR. 

A monitoring programme with this telescope on extragalac- 
tic sources visible by both the NRT and H.E.S.S. is in place 
since 2001. For the campaign described here it consisted of a 
measurement at 1 1 cm every two or three days. The average 
flux for the 8 measurements in October and November 2003 
was 0.30 ± 0.01 Jy with possible marginal variability. 

3. Results 

3.1. Light curves 

The October and November 2003 light curves of all the 
H.E.S.S., RXTE and ROTSE observations are shown in 
Figure [2] The H.E.S.S. light curve is binned in run-length 
times averaging 28 minutes each. The flux is in units of 
10 -11 photons cm^V above 300 GeV, derived using the av- 
erage photon index 3.37 obtained in section l3~3l Spectra could 
not be derived on a run-by-run basis due to the weak signal. 
As for the observations reported in AH04, the overall light 
curve is inconsistent with a constant flux. A x 2 fit of the data 
to a constant yields a 3.4 x 10~ 10 x 2 probability. The intra- 
night VHE flux on MJD52936 (Fig. [3} exhibits an increase of 
a factor of 1.9 + 0.6 in 0.1 1 d. On MJD52932 the peak-to-peak 
flux shows an increase of a factor of 2.5 + 0.9 within 0.09 d. 
These timescales are longer than the 30 minute doubling time 
reported in AH04. For these two extreme cases of VHE vari- 
ability observed during this campaign only the second had a 
limited RXTE coverage. 

The 2-10 keV X-ray flux in this campaign ranges from 
-FViokeV = 2.0x 10~ n ergcm^s -1 to 4.4x 10~ n erg cirT 2 s" 1 . 
The maximum is lower than the 20 November 1997 measure- 
ment of 2.3 x 10" 10 erg cm -2 s _1 dVestrand~& Sreekumar 1999 1 
indicating that the X-ray state seen here is not exceptionally 
high. The minimum seen here is consistent with historically 
low fluxes jZhang, Treves et al. 20 02 1. The intra-night variabil- 
ity is also obvious here, but no flare was completely resolved. 

A 60% flux variability in f val - * 1.5 ks on MJD52936 is 
the best marked transient episode in the observations reported 
here (bottom panel b) in Figure |3]for which the H.E.S.S. ob- 
servations were made at the end of the transit inducing a large 
associated error on the flux estimation due to the high zenith 
angle of the source. This timescale is comparable to those re- 
ported by Gai dos et al. (1996)| where doubling times as short as 
15 min from Mkn 421 were observed in the VHE band. 

This flare is the fastest rise seen in this object to 
date since BeppoSAX saw a 5 x 10 4 s rise timescale 
(Zha ng, Treves et al. 200 2 1 and Katao ka et al. (2 000) observed 
a doubling timescale of 3 x 10 4 s with the ASCA satel- 
lite. So far the fastest rise in this type of object was 
observed in Mkn 501 with a 60% increase in less than 
200 s JCatanese & Sambruna 20001 though |Xue"& Cui (2004)] 
claims that this flare is likely to be an artifact. 
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Fig. 2. Top: (a) October 2003 Light curve from H.E.S.S. binned in run lengths averaging 28 minutes each, in units of integral flux 
above 300 GeV. (b) Light curve expressed as flux in the 2-10 keV band for the RXTE observations, in 400 s bins, (c) Light curve 
derived from the ROTSE optical data. Bottom: Same as above but for the November 2003 data. Note that the RXTE observations 
were much shorter than in the previous month. 
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The optical emission of PKS 2155-304 is dominated by 
the nucleus which outshines the host galaxy by a factor of « 4 
given in KFS98. The observed variability amplitude is there- 
fore not biased by the constant emission of the host galaxy and 
mainly due to the activity of the nucleus. The peak-to-peak am- 
plitude of variability is moderate compared to the variability 
amplitude at shorter wavelengths and typically 0.1 mas peak-to- 
peak. The object has been monitored over longer time-scales 
with the ROTSE-IIIc telescopes showing variations with am- 
plitudes close to l ma s. 
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Fig. 3. (a) The VHE light curve derived run by run on 
MJD52936. The horizontal error bars are the length of the run 
from which the flux is derived. The dashed line is the result of a 
X 1 fit of the data to a constant, which yields a x 2 of 15 for 7 de- 
grees of freedom, corresponding to a 0.04 x 1 probability which 
is evidence for the variability seen here, (b): The 2-10 keV X- 
ray flux {open circles, right scale) and counting rate normalized 
to 1 PCU (lines, left scale) showing the fast transient. 
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Fig. 4. Correlation plots between different wavebands. No 
clear correlation is found between any of the simultaneous 
measurements, (a) Correlation plot for the 23 X-ray data 
segments that overlapped exactly with a H.E.S.S. observa- 
tion, (b) Simultaneous H.E.S.S. and ROTSE observations, (c) 
Simultaneous RXTE and ROTSE observations. The optical 
data were binned to the 400 s long RXTE segments, which are 
overlapping on 6 different days. 



3.2. Correlation analysis 

In order to quantitatively look for correlated variability be- 
tween the VHE, X-ray and optical bands, the measured fluxes 
are plotted against each other in Figure @] for all the observa- 



tions carried out during this campaign. For correlated VHE/X- 
ray variability, the RXTE analysis was slightly modified: only 
observation segments that happen exactly within a H.E.S.S. run 
are reduced and analyzed (using only PCU2, PCU3 and combi- 
nations thereof). This provides 23 simultaneous data segments 
for the whole campaign for which the fluxes are represented in 
Figure^ There is no obvious correlation for those observations 
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Fig. 5. Time-averaged spectrum derived from the October and November 2003 H.E.S.S. data along with a fit to a powerlaw. 



(correlation factor r = 0.37 ±0.1 3). For correlated optical/VHE 
variability, ROTSE observations that happen within a H.E.S.S. 
run are averaged and their errors summed quadratic ally. No 
correlation (r = 0.24 ± 0.27) is found for these observations. 
Also no correlation was found between the optical and X-ray 
band (r = -0.02 ± 0.05). 

3.3. Spectra 

The October and November 2003 H.E.S.S. data were all com- 
bined for the spectrum shown in Figure |5] The best-fitting 
power law (x 2 = 7 for 8 degrees of freedom) is given by: 



§ = ,2,3 ± 0,7>xl0-(4) 



-3.37±0.07±0.10 



cm 2 s 



TeV" 



(1) 



whichis comparable to F = 3.32 ±0.06 and / = (1.96±0.12)x 
10~ 12 cirT 2 s _1 TeV -1 previously reported in AH04. 

The result of the broken power law fit for the combined 
RXTE PCU2 and PCU3 spectrum is shown in Figure |6] It 
yields an unabsorbed flux in the 2-10 keV band of .F2-10 kev = 
(2.66 ± 0.04) x 10" 11 erg ctrrV 1 Of 2 = 41, 31 degrees of free- 
dom), a lower index of Ti = 2.81 ± 0.05, a break energy of 
E b = 4.9 ± 0.8 keV and a higher index of T H = 2.95 ± 0.04. 
A single power law fit to the same data yields F2-10 kev = 
(2.69 ± 0.03) x lO^'ergcirrV 1 (y 2 = 51, 34 degrees of 



freedom) and an index of F = 2.88 ± 0.13, a poorer fit than 
the broken power law, but the index still provides information 
that can be used for comparison with historical measurements. 
Indeed, the derived index is close to those measured by the 
BeppoS AX satellite ( Giommi, Fiore, Guaina zzi et al. 19 98 1, 
GINGA jSembay et al. 1993j and well w ithin the range 
observed bv EXOSAT JTreves et al. 1989b . The statistics 
above lOkeV in our RXTE observations are too poor to 
check the existence of a possible hard tail above 20keV 
( |Giommi, Fiore, Guainazzi et al. 1 998 ) which might be the sig- 
nature of the onset of a high-energy component. 

In order to look for flux dependent spectral variability, the 
RXTE data subset used in Figure [5] is divided into two en- 
ergy bands, the PHA channels 0-9 (soft band) and 10-27 (hard 
band), corresponding to approximately 1-4 keV and 4-1 1 keV, 
respectively. A hardness ratio (HR), shown in Figure is the 
ratio of the counting rate in the hard band over the soft band. 
There is a clear correlation of the HR with the rate, peaking 
when the rate is highest. The correlation factor between the 
rate and the HR is r — 0.76 ± 0.12. Even though the variabil- 
ity timescale here is much smaller, this behavior is compatible 
with the hardening reported in Chiapp etti et al. (1999)| 

3.4. EBL corrected spectrum 

For objects at non-negligible redshifts, the large, energy- 
dependent opacities can cause the emitted spectrum to 
be greatly modified both in shape and intensity (see 
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Fig. 6. X-ray spectrum derived from the summed October and November 2003 data, using PCU2 and PCU3, fitted by a broken 
power law. The top panel shows the data and folded model, the bottom panel shows the residuals between the data and the model. 




2.5 3 3.5 4 4.5 
Rate (Counts/s/PCU) 

Fig. 7. Plot of the hardness ratio HR versus the counting rate 
normalized to 1 PCU. 



e.g.lStecker et al. 1992llBiller 1995l|Coppi & Aharonian 1999| 
Vassiliev 2000 1. Unfortunately, at present the knowledge of the 
EBL still has large uncertainties, for both direct measurements 
and models, as summarized in Prima cFet al. (200"T)| In order 
to estimate the intrinsic VHE spectrum, and thus to locate the 
Inverse Compton (IC) peak of the blazar's SED, we have used 
three EBL models (Fig. [HJ as representatives of three different 
flux levels for the stellar peak component. This is the EBL en- 
ergy range which mostly affects the H.E.S.S. spectrum: with 
data up to 3-3.5 TeV, the peak of the y — y cross section is 
reached for soft photons with wavelengths < 5 fim. 



The three models used here are (in order from higher 
to lower fluxes): the phenomenological shape used in 
Aharonian et al. (2003) which is based on the original 
Primacket al. (2001) calculation but smoothed and scaled up 
to match the data points below 1 fim and at 2-3.5 /im (here- 
after Phigh); the original Primack et al. (2001) calculation 
for a Salpeter initial mass function (hereafter PrimackOl); 
and the new 2004 calculation ( Prim ack et al. 20 04- hereafter 
Primack04), which takes advantage of the recent improvements 
in the knowledge of the cosmological parameters and of the lo- 
cal luminosity function of galaxies. 

The opacities are calculated from the EBL SED shapes tak- 
ing into account only the cosmology (Hq = 70km/s/Mpc, 
£^Mat = 0.3 and Qa = 0.7). To treat all three shapes simi- 
5 larly, no evolution has been introduced at this point. This cor- 
responds to a "maximum absorption" hypothesis (i.e., for in- 
creasing z, constant instead of decreasing EBL comoving en- 
ergy density). But at these redshifts (~ 0.1) and, for example, 
assuming the evolution given in Prim ack etal. (200 1)| the dif- 
ference is still small (variation in the photon index AT < 0.1 
in the range 0.3-1 TeV), and negligible compared to the differ- 
ences between models. 

The resulting absorption-corrected spectra are shown in 
Figure [9] together with the observed spectrum. The intrinsic 
spectra are all well fitted by a single power-law model, with 
a hard spectrum for Phigh (F * 1.5), and soft spectra for 
PrimackOl and Primack04 models (r ss 2.3 and T as 2.8, re- 
spectively). This effect is directly related to the different flux 
levels of the stellar peak component, which imprint a differ- 
ent amount of softening onto the original spectrum. This di- 
rect link thus yields two simple scenarios for the location of 
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Fig. 9. PKS 2155-304 absorption-corrected spectra, along with the observed spectrum, for each EBL model considered and 
in a vF v plot (i.e. a zoom in the blazar's SED). Open dots: H.E.S.S. observed spectrum. Filled dots: reconstructed intrinsic 
spectrum. All spectra from "Primack-type" models are compatible with a power-law shape, with different slopes: hard (Phigh) 
and soft (PrimackOl, Primack04), locating the intrinsic SED peak for the high component above 1-2 TeV or below 200 GeV, 
respectively. 



the blazar's high energy peak, with the dividing line repre- 
sented by the EBL flux which gives F = 2.0 (1.3 times the 
PrimackOl model). Models with stellar peak fluxes above this 
(such as Phigh, and generally all those in agreement with the 
direct estimates of the fluxes between 2 and 3.5 /mi) imply a 
hard intrinsic spectrum (r < 2), and thus an IC peak above 1— 
2 TeV. EBL models with lower fluxes (such as the PrimackOl 
and Primack04) imply instead a soft spectrum (T > 2), locating 
the IC peak below the observed energy range (< 200 GeV). In 
the following, we will discuss both scenarios for the SED mod- 
elling, using the Phigh and Primack2004 curves as the two ends 
of the possible range of values for the "Primack-type" shape 
(i.e., between the claimed EBL direct measurements at few mi- 
crons and the lower limits from galaxy counts). 

4. SED Modelling 

The broadband spectral morphology of PKS 2155-304 is typi- 
cal of the BL Lac type, with a double-humped structure in vF v 
representation, exhibiting a low-energy and a high-energy com- 
ponent. Its broadband emission is usually attributed to emis- 
sion from a beamed relativistic jet, oriented in a direction 
close to the line of sight ( Begelman, Blandfor d & Rees 1984| 
Blan dford & Konigl 1979) . The spectral energy distribution in 
units of power per logarithmic bandwidth E 2 dN/dE versus en- 



ergy E is shown in Figure ^3 The EGRET measurements, 
between the H.E.S.S. and RXTE points, are from the third 
EGRET catalog (Hartman 1999 1 and from a very high y-ray 
state described in |Vestrand, Stacy & Sreekumar (1995)| There 
is a difference in spectral states, since in the former case the 
power law photon index is 1.71+ 0.24 whereas it is 2.34 ± 0.26 
in the latter which most likely consists of a mix of low and 
high activity state observations. The historical EGRET spectra 
are therefore unlikely to represent the state of PKS 2155-304 
during the campaign presented here and are not used to put 
stringent constraints on the modelling. Considering the archival 
data and the steep X-ray spectrum in Figure [lO] the peak 
of the low-energy component occurs in the 2-2000 eV range. 
The archival BeppoSAX data from a high state analyzed by 
Chiappet ti et al. (19 99) and represented here above our RXTE 
data, show a peak at » 0. 1 keV. The absorbed VHE peak loca- 
tion is clearly below 300 GeV, with its exact location depend- 
ing on the spectrum in the EGRET range. 

Whereas the current models seem to agree that the 
low-energy component is dominated by synchrotron ra- 
diation coming from nonthermal electrons emitted in 
collimated jets, the high-energy emission is assumed 
to be either inverse Compton scattering off the syn- 
chrotron photons (Synchrotron Self-Compton, SSC, see 
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Fig. 8. Spectral energy distribution of the Extragalactic 
Background Light (EBL). It is thought to be characterized 
by two distinct bumps, around 1-2 and 100-200 /vm, pro- 
duced by the red-shifted stellar light and re-radiated ther- 
mal dust emission, respectively. The compilation of mea- 
surements have been taken from |Hauser & Dwek (200 1)| and 
|Aharonian et al. (2003)| Lower limits are from HST and 
ISO source counts JMadau & Pozzetti 20001 [Elbaz et al. 20021 
IGardner et al. 20001 Above 400 jj.m, the contribution of the 
CMB starts to dominate. The three continuous lines represent 
the "Primack-type" EBL shapes: from higher to lower fluxes 
(at 1 yum), Phigh, PrimackOl and Primack04 models (see text). 
The steep line starting at 300 /urn is the onset of the CMBR. 



e.g. Marasc hietal. 19991 IBicknell & Wagner 2002) or by 
external photons (see e.g. ISikoraetal. 1994i . This kind of 
leptonic model will be discussed in section 14.21 A hadronic 
origin of the VHE emission using the Synchrotron-Proton 
Blazar (SPB) model with a dominating proton synchrotron 
component at high energies in a proton-electron plasma is 
also able to produce a double humped SED and is discussed 
in section l4~3l The lack of correlation between the RXTE and 
H.E.S.S. fluxes (and possibly also between the optical and 
the VHE emission) within the small variability range may 
indicate a different spatial origin, or a different underlying 
particle distribution. In the proton synchrotron model a lack 
of correlated variability between y-ray emission and the low 
energy electron synchrotron component could arise if the 
electrons and protons are not co-accelerated. 

The high-energy component above » 100 GeV is attenu- 
ated by interactions with the EBL and is a lower limit for the 
intrinsic spectrum. The energy budget in X-rays and VHE y- 
rays is comparable, though the maximum output at the peak 
energy in the high-energy component is likely to be lower than 
that in the lower-energy component. Interpolating between the 
high-state EGRET archival data and VHE data would lead to a 
maximum located above 10 GeV, which is surprising since the 
observations reported here indicate a low state. Extrapolating 
the EGRET catalogue spectrum to VHE energies with a power 



law falls below the H.E.S.S. data and therefore requires two 
inflexion points in the SED. Simultaneous observations in the 
MeV-GeV range with the upcoming satellite GLAST will be 
crucial to constrain the high-energy component shape. 

4.1. Doppier boosting and synchrotron/Compton 
derived parameters 

The electromagnetic emission in blazars is very likely to 
be Doppler-boosted (or beamed) toward the observer. In 
the radio regime, the evidence for Doppier boosting in 
PKS 2155-304 comes from superluminal expansions observed 
with VLBI (Piner & Edwards 2004 1. Relativistic beaming is 
also required in order to avoid absorption of GeV photons by 
X-ray photons via the e + /e~ pair-production process (see e.g. 
|Maraschi, Ghisellini & Celotti 1992} . It is thus possible to use 
the y-ray variability to establish a limit for the Doppier factor 6, 
with 5 defined in the standard way as |T(1 -yScos ff)]~ l , where 
T is the bulk Lorentz factor of the plasma in the jet, /? = v/c, 
and 9 is the angle to the line of sight. 

Following Dondi & Ghisellini (1995) the size of the y-ray 
emission zone R is derived from the time-scale of variability 
? V ar (supposing the timescale of intrinsic variability is negligible 
compared to the light crossing time) by R < cf var 5/(l + z). In 
this case (assuming a time scale f var as 2 ks from section^) the 
size of the emission region is 



R6 1 < 5 x 10° cm. 



(2) 



Since for y-rays in the 300 Ge V-3 Te V range the opacity 
T yy cannot significantly exceed 1, independently of the emis- 
sion mechanisms, another constraint on the minimum value 
of the Doppier factor can be derived by estimating how much 
Doppier boosting is necessary for photons with observed en- 
ergy E y to escape from a source with radius R and a flux den- 
sity F(E,), where E t = (m e c 2 ) 2 6 2 /Ey(l + zj 1 . At this point 
one can follow [Mattoxet al. (1993)| and especially Eq. 3.7 in 
Pondi & Ghisellmi^995)| writing that the opacity is 



(Tj 



i £ 



l 



5 he R 6 3 (l +z) 



F(E t ). 



(3) 



Imposing that r yr < 1 for photons with observed energy E y = 
1 TeV yields a lower limit on 6 for given R that can be derived 
numerically from the observed SED. For PKS 2155-304 this 
yields 

R 1 (T 6 - 4 < 5 .6 x 1 0~ 24 cm" 1 . (4) 

The minimum allowable boost 6 comes from combining Eq.|2] 
and Eq.|4]and yields 

S>19 (5) 

which is higher than the limit obtained in a similar way in 
|Tavecchio et al. (199 8) but at the lower end of the range usu- 
ally obtained from SSC modelling as in |Kataoka et al~ 2000) 
It is useful to stress that this constraint is valid under the usual 
assumption that the region emitting the SED optical flux (i.e., 
the target photons) is cospatial with (or at least embeds) the 
high energy emitting region. The models used later in this pa- 
per are therefore essentially single zone models. 
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If the observed X-rays are synchrotron radiation from non- 
thermal electrons then the mean observed energy < E > of an 
electron with Lorentz factor y e is given by 

6 2\%ylqB 

<E>* - 

1 + z 15 V3m e 

Using < E >= lOkeV here yields 

B6y 2 e = 1.1 x 10 12 G. (6) 

In the SSC scenario the same electrons can Comptonize ambi- 
ent photons up to the VHE regime, thus 

y e m e c 2 > 3 TeV. (7) 

1 +z 

Combining Eq.|S]and Eq.Qyields 

— %-) (i + zr 2 . (8) 

m e c z ) 

A numerical application yields 

B8~ l < 0.03 G (9) 

and hence the Lorentz factor is constrained to y e > 1.3 x 10 6 
which is higher than that calculated as above for Mkn 421 
bv lTakahashietal. 19961 who derived y e > 5 x 10 5 (and B * 
0.2 G) but lower than y e ~ 10 7 found in a similar way in 
1ES 1959+650 Kiiebels et al. 20021 . 

The X-ray data presented above imply that the X-ray spec- 
trum of PKS 2155-304 hardens as the source brightens. This 
is often measured in BL Lac objects; a hardening of the spec- 
trum when flares occur, and a blueward shift of the peak of 
the synchrotron emission v sync (and presumably higher energy 
inverse-Compton emission) by factors that can be as large as 
100 were measured in the cases of Mkn 501 (Pian et al. 1998t . 
1ES 1426+428 and PKS 0548- 322 dCostamante et al. 20011 . 
In the case of PKS 2005-489 JPerlman et al. 19991 . a more 
moderate shift of a factor of 3 or less of the synchrotron emis- 
sion was found. The archival data suggest that v sync lies in the 
UV band for PKS 2155-304, but no data were taken simulta- 
neously in this campaign at that wavelength. 

The lack of indication for correlated X-ray/VHE variability 
does not imply that PKS 2155-304 behaves differently from 
VHE sources such as Mkn 421, for which VHE/X-ray corre- 
lation has been established on a much higher variability basis 
(see e.g. lCui et al. 2 004 1 with dynamical ranges of 30 in both 
energy bands. Limiting those observations to the same dynam- 
ical range as observed here would not allow any claim for cor- 
relation. Future observations of PKS 2155-304 with a higher 
variability amplitude would bring more insight into this. 

4.2. Leptonic interpretation 

Interpretation with a single zone SSC model of the SED 
of PKS 2155-304 has already been proposed in the litera- 
ture using two different assumptions. In Katao ka et al. (2 000) 
the low energy tail of the SSC model is used to account 
for the low-energy component up to the optical in the SED. 



Parameter 


Model 1 


Model 2 


Model 2 Phigh 


fl blob (10 15 cm) 


6.5 


1.5 


1.5 


B(G) 


0.15 


0.25 


0.02 


S b 


25 


25 


50 


K(cm- 3 ) 


160 


2.0 x 10 3 


240 


7 brea k(10 3 ) 


7.5 


100 


300 


7max(10 5 ) 


3.8 


9.0 


5.0 


»1 


1.4 


1.7 


1.6 




3.2 


4.65 


4.5 


fl Jct (10 6 cm) 


1.0 


1.0 


1.0 


^jet(pc) 


55 


55 


55 


<5jet 


2 


2 


2 


ifjet(cnr 3 ) 


40 


40 


40 


Sje,(G) 


0.04 


0.04 


0.04 


7breakjet(10 3 ) 


2.5 


45 


45 



Table 1. Parameters for the SED fit with the assumption that 
the optical and X-ray emission are part of the jet synchrotron 
emission (Model 1 in the table, dotted line in Fig. 1 1 Oi or that 
the optical emission emanates from the VLBI emission zone 
(Model 2 in the table, dashed lines in the same Figure). The 
parameters are described in |Katarzynski et al. (200 1)| 



That component is decomposed into two sub-components by 
Chiappet ti et al. (1999)| where the radio to optical emission 
has another origin than the X-rays, which are assumed to 
come from the jet. These two different interpretations are 
used here in the context of the leptonic model described in 
Katarzyhski et al. (2001) which has already been applied to 
Mkn 501 and Mkn 421 l |Katarzynski et al. 2003) . To constrain 
this model, only the simultaneous data are used, since the 
archival data reported in the SED of Figure ^3 &g likely to 
not represent the state of this source (note for example the dif- 
ference in optical flux and the ROTSE measurement). 

When using the Primack04 absorption, the model used here 
can reproduce the X-ray through VHE part of the SED, but the 
H.E.S.S. spectrum constrains it such that the radio measure- 
ment can not be included in the synchrotron bump predicted by 
the single-zone model. As for Mkn 421 and Mkn 501, adding 
a more extended component than the VHE emitting zone can 
provide an explanation for this. The origin is probably the 
compact VLBI core which has a radio core to lobe ratio of 
* 1 ULaurent-Muehleisen et al. 19931 IPiner & EdwardT2 004 1 
and a typical size of 10 18 cm, more than two orders of mag- 
nitude larger than the VHE emitting zone. This VLBI fea- 
ture dominates the spectrum at low energy and is included 
in the SED modelling here. An uncertainty remains which 
is the high frequency cutoff of this VLBI component. The 
host galaxy contribution to the optical flux is estimated to be 
w 10~ n ergcm~ 2 s _1 , deduced from the magnitudes given by 
KFS98 and assuming a low-redshift solar metallicity elliptical 
galaxy of age equal to 13Gyr {R - H = 2.4), corresponding 
to a mass of 5 x 10 n M o ( |Fioc & Rocca-Volmerange 1997) . So 
even at the measured low activity state of PKS 2155-304 the 
host galaxy is not contributing much in the optical range. 

The ROTSE measurement can be ascribed here to either the 
high-energy tail of the VLBI component or to the synchrotron 
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Fig. 10. Spectral energy distribution of PKS2 1 55 -304. Only simultaneous measurements are labeled. Non contemporaneous data 
are in grey symbols. The H.E.S.S. spectrum is derived from October and November 2003 data (filled circles) as is the RXTE 
spectrum. The NRT radio point (filled square) is the average value for the observations carried out during this period. The two 
triangles are the highest and lowest ROTSE measurements for the Oct-Nov observations. Archival SAX data represent the high 
state observed in 1997 from Chiappetti et al. (1999) Archival EGRET data are from the third EGRET catalogue (shaded bowtie), 
Hartman 1999 1 and from a very high y-ray state described in Vestrand, Stacy & Sreekumar (1995) (open bowtie). 1995 U B V 
R I data are from Paltani ^t al. (1997)| where the mean flux and flux deviation is used here. IUE data from |Urry et al. (1993)| at 
2800 A and 1400 A are included here with an error bar corresponding to the spread of the measured fluxes. Other data are NED 
archival data. The solid line is the hadronic model described in section l4~3l The dotted and dashed lines are the same leptonic 
model with different assumptions described in section l4~2l the dotted line assumes a common origin for the optical and X-ray 
synchrotron emission, while the dashed line is the case where the optical emission emanates from the VLBI core; for both cases 
the lowest frequency hump is the predicted VLBI core emission. All the VHE emission in the models are absorbed according to 
the Primack04 model discussed in section l3~4l 



part of the SSC model. Assuming a common origin for the X- 
rays and the optical emission, and using a variability time scale 
of 0. 1 d to constrain the emitting zone, the model tends to pre- 
dict a high IC flux, as shown in Figure [TO] However, the lack 
of correlation between the X-rays and the optical emission in 
our measurements - also suggested by Domin ici et al. (20 04) 
based on less sensitive RXTE/ASM measurements - indicates 
that the optical emission could originate from the VLBI com- 
ponent, which is modelled by a slight increase in the maximal 
Lorentz factor of the emitting electrons. This in turn lessens the 
constraint on the simultaneous SSC fit of the X-ray and VHE 
part and allows a better fit of the VHE spectrum for smaller 



sizes of the emitting zone. Detailed parameters of the two hy- 
potheses are given in Tabled 

If the absorption correction is well described by the Phigh 
model, the slope of H.E.S.S. data at high energy implies that 
the peak of the TeV emission bump is located above 4TeV 
(or 10 27 Hz). Such a high frequency peak emission imposes a 
strong constraint for the single-zone SSC scenario, especially 
since the peak of the synchrotron bump has to remain below 
1 keV (or 10 17 5 Hz) as required by the slope of the RXTE data. 
High values of both the jet Doppler factor and the maximal 
Lorentz factor of radiating particles are required to reach the 
necessary energy for the IC bump, that is y max c> > 8 x 10 6 . 
On the other hand, to keep the synchrotron peak below 1 keV 
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imposes an upper limit to the magnetic field. Within these con- 
strains, the best fit we obtain is shown in Figure ^2 We can 
note that none of the high energy tails are well accounted for. 
The set of parameters for the best fit is given in Tabled This fit 
marginally reproduces the observed X-ray and y-ray data, but 
is not as satisfactory than that obtained with the Primack04 ab- 
sorption correction, and in any case it is impossible to take into 
account the ROTSE optical point. The main changes in param- 
eters between the two fits consist in enhancing the boosting, 
which then becomes quite extreme, while reducing the density 
and magnetic field for the Phigh absorption correction. 

The constraints on 5 derived here from either simple opac- 
ity arguments or from the one-zone model parametrisation of 
the SED are in the range of Doppler factors usually derived 
with such assumptions or models for other VHE emitters. As 
pointed out by |Chiaberge et al. (2000)| such high values are 
however at odds with attempts to unify the BL Lac popula- 
tion with the family of FR I sources ( |Urry & Padovani T9 95 1, 
the latter being possibly an unbeamed since off-axis viewed 
case of the former. The same authors suggest that mod- 
els where velocity structures in the jet, such as the "spine- 
sheath" model (see e.g. |Sol, Pelletier & Asseo 19891 or the 
decelerating flow model flp eorganopoulos & Kazanas 2003 
|Geor ganopou los & Kazanas 2004) allow lower bu lk Lorentz 
factors. Another option from Pelleti er et al. 20 04 is to cope 
with the pair creation catastrophe implied by smaller Doppler 
factors in their "two-flow" solution. Comparing the SED with 
such models, which make the BL Lac - FR-I connection more 
plausible, is an interesting task but beyond the scope of this 
paper. 

4.3. Hadronic models 

Generally, the leptonic models constitute the preferred con- 
cept for TeV blazars, essentially because of two attrac- 
tive features: (i) the capability of the (relatively) well un- 
derstood shocks to accelerate electrons to TeV energies 
( |Sikora & Madejski 2001| IPelletier 200H and (ii) the effec- 
tive conversion of the kinetic energy of these relativistic 
electrons into the X-ray and VHE y-ray emission compo- 
nents through the synchrotron and inverse Compton radia- 
tion channels. The hadronic models are generally lack these 
virtues. They assume that the observed y-ray emission is ini- 
tiated by accelerated protons interacting with ambient matter 
JBednarek 1993IIDar & Laor 1997IIPohl & Schlickeiser 2000L 
photon fields (PIC model, Mannheim 2000 1, magnetic fields 
JAharonian et al. 2002b or both JMiicke & Protheroe 20"00Ti . 

The models of TeV blazars involving interactions of pro- 
tons with photon and B-fields require particle acceleration to 
extreme energies exceeding 10 19 eV which is possible if the 
acceleration time is close to f acc = r g /c (r g is the gyro-radius). 
This corresponds (independent of a specific acceleration mech- 
anism) to the maximum (theoretically possible) acceleration 
rates ( Aharo nian et al. 2 002 1 which can only be achieved by 
the conventional diffusive shock acceleration in the Bohm dif- 
fusion regime. 



On the other hand, the condition of high efficiency of radia- 
tive cooling of accelerated particles requires extreme parame- 
ters characterizing the sub-parsec jets and their environments, 
in particular very high densities of the thermal plasma, radia- 
tion and/or B-fields. In particular, the proton-synchrotron mod- 
els of TeV blazars require highly magnetized (B » 10G) con- 
densations of y-ray emitting clouds containing Extremely High 
Energy (EHE) protons, where the magnetic pressure dominates 
over the pressure of relativistic protons (Aharonian 2000 1. 

Below we use the hadronic SPB model 
(Mu cke & Protheroe 20001 IMiicke & Protheroe 200 fl to 
model the average spectral energy distribution (SED) of 
PKS 2155-304 in October-November 2003. A detailed 
description of the model, and its implementation as a (time- 
independent) Monte-Carlo/numerical code, has been given in 
IMiicke, Protheroe et al. (2003)| and |Reimer et al. (2004)| 

Considering the rather quiet activity state of 
PKS 2155-304 in Oct-Nov. 2003, we use the 3EG cata- 
log spectrum, since it is the best determined EGRET spectrum 
from this source to date, as an upper limit for modelling 
purposes. 

Flux variability provides an upper limit for the size of the 
emission region. To allow for a comparative study between lep- 
tonic and hadronic models we fix here the comoving emission 
region to R ~ cf val -<S = 5 x 10 13 5cm deduced from the X-ray 
variability. We assume that the optical through X-ray emission 
and the y-ray output stem from the same region of size R. 

A reasonable model representation for the simultaneous 
data assuming a Primack04 model for the VHE absorption is 
found for the following parameters: magnetic field B = 40 G, 
Doppler factor 6 = 20, injection electron spectral index a e = 
1 .6, assumed to be identical to the injection proton spectral in- 
dex Cp, maximum proton energy of order y p , max ~ 4 x 10 9 , 
e/p-ratio of 0. 15 and a near-equipartition proton energy density 
of u p = 27 erg crrT 3 . The required total jet power is of the order 
Ljet ~ 1.6 X 10 45 erg s _1 . When using the Phigh EBL model, a 
reasonable representation of the data may be achieved by in- 
creasing the maximum injected proton energy to y Pjmax = 10 10 
and simultaneously increasing the e/p-ratio to 0.24, while all 
other parameters remain unchanged. Note that here the max- 
imum proton gyro-radius approaches the size of the emission 
region. Alternatively, a doubling of the magnetic field to 80 G 
together with an increase of y Pjlnax to 8 x 10 9 and a e/p-ratio of 
unity (leading to Kg « 50m p ) represents the SED-data equally 
well. In conclusion, none of the "Primack-type" EBL mod- 
els can explicitly be ruled out in the framework of the SPB- 
model by the H.E.S.S. data presented here. In all cases, proton 
synchrotron emission dominates the (sub-)TeV radiative out- 
put. Depending on the Doppler factor, part of the proton syn- 
chrotron radiation produced may be reprocessed to lower en- 
ergies. Contributions from the muon and pion cascades are al- 
ways lower than the proton synchrotron component. The low 
energy component is dominated by synchrotron radiation from 
the primary electrons, with a negligible contribution of syn- 
chrotron radiation from secondary electrons (produced by the 
p- and /^-synchrotron cascade). 

On the other hand, the synchrotron radiation of secondary 
electrons resulting from interactions of VHE y-rays with ex- 
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Fig. 11. Estimations of the intrinsic PKS 2155-304 VHE spec- 
trum, corrected for either the Phigh or Primack04 models, 
along with the associated intrinsic (i.e. before absorption) SSC 
model (dashed lines) and SPB model (solid lines). 



ternal low-energy photons with a modest yy — > e + /e~ opacity 
(r yr < 1) may lead to significant X-ray emission with a lumi- 
nosity comparable to the luminosity of the primary VHE emis- 
sion (Aharonian 2000). 

Models involving meson production inevitably predict neu- 
trino emission due to the decay of charged mesons. The SPB- 
model for PKS 2155-304 explains the high energy emission 
dominantly as proton synchrotron radiation, making the neu- 
trino flux completely negligible. 

5. Conclusions 

This paper reports multi-wavelength observations of the BL 
Lac object PKS 2155-304 in 2003. Although the source ap- 
peared variable in the VHE, X-rays and optical bands, the latter 
two indicate that PKS 2155-304 was in a state close to histor- 
ically low measurements, even though it was easily detectable 
by H.E.S.S. in all nights of observations since the beginning 
of the detector operation (see AH04 for the observation his- 
tory up to August 2003). An extreme case of VHE variability 
shows a peak-to-peak increase of a factor of 2.5 + 0.9 in 0.09 d. 
Variability on the timescale of a few ks in the 2-10 keV band 
and of the order of 0. 1 d for energies > 300 GeV were observed 
by RXTE and H.E.S.S. The X-ray data show a correlation of 
the flux with the spectrum, which becomes harder when the 
source is brighter. At the level of the simultaneously observed 
modest variability, no correlation between the VHE y-rays, the 
X-rays and optical was seen. Observations with greater vari- 
ability and better coverage are needed before it can be asserted 
that the VHE/X-ray pattern in PKS 2155-304 is different from 
other known VHE-emitting AGN. Since the source was in a 



low emission state in both the optical and X-rays compared 
to archival measurements, this lack of correlation has yet to 
be established for a higher emission state. Simultaneous obser- 
vations in the X-rays/optical band and VHE y-rays had never 
previously been performed on this scale. Its continual VHE 
detection makes PKS 2155-304 unique in the TeV BL Lac 
category, and probably indicates that H.E.S.S. has achieved 
the sensitivity level where it can detect the quiescent state of 
PKS 2155-304 at any time. A time-averaged energy spectrum 
is determined for the 2 observation periods and fits to a power 
law (F = 3.37 + 0.07) in the VHE y-rays, and to a broken power 
law (T L = 2.81 ± 0.05, E b = 4.9 ± 0.8 keV, Y H = 2.95 + 0.04) 
in the X-rays. 

A comparison of the intrinsic spectrum with predictions 
from existing one-zone leptonic and one-zone hadronic mod- 
els is attempted to give a plausible estimation of underlying 
physical parameters. The values of the parameters are in line 
with those inferred for other VHE-emitting blazars. In these 
models the VHE emission is attenuated according to two dif- 
ferent EBL levels. This changes mainly the Doppler boosting 
in the leptonic model, but the high level EBL decreases the 
agreement significantly. In the hadronic model, the maximum 
injected proton energy can be changed to accomodate different 
EBL levels and can therefore not significantly constrain any of 
the EBL models used here. 
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